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ABSTRACT 


_ S 

It  is  shown,  with  the  aid  of  two  examples,  that  the  calculus  of 
variations  may  be  used  to  an  advantage  over  the  method  of  matched 
asymptotic  expansions  in  solving  singular  boundary-layer  problems  of 
applied  mechanics. 


INTRODUCTION 


There  are  only  few  methods  available  to  solve  singular  boundary-layer 
problems  in  applied  mechanics  [1]  to  [2].  One  of  them  is  the  method  of 
matched  asymptotic  expansions. 

The  problem  usually  consists  in  finding  a  solution  to  certain  dif¬ 
ferential  equations  (subject  to  appropriate  boundary  conditions),  in  which 
the  highest  derivative  is  multiplied  by  a  small  parameter.  One  then  tries 
to  use  a  regular  perturbation  scheme  to  solve  the  problem  and  finds  sev¬ 
eral  solutions,  none  of  them  being  valid  over  the  entire  region  of  interest. 
However,  these  solutions  usually  contain  some  unknown  constants  which  are 
determined  by  matching  those  solutions  with  each  other  over  their  common 
domain  of  validity. 

In  what  follows,  we  will  show  with  the  aid  of  two  boundary  value  prob¬ 
lems,  that,  instead  of  using  the  matched  asymptotic  expansion  method  to 
find  those  unknown  constants,  the  calculus  of  variation  may  be  used.  This 
means  that  the  calculating  procedure  is  thus  reduced  to  a  systematic 
series  of  steps  not  requiring  any  particular  skill  or  experience.  Of  course 
the  procedure  to  be  described  may  be  used  only  in  those  problems  where  a 
variational  formulation  may  be  made  effected,  such  as  in  potential  flow 
problems  in  fluid  dynamics  and  in  conservative  problems  of  solid-body 
mechanics . 

The  two  examples  considered  here  are:  1)  Cylindrical  deflection  of  a 
slightly  stiff  membrane,  and  2)  Torsion  of  a  shaft  with  a  keyway.  The 
first  problem  is  described  by  an  ordinary  differential  equation,  whereas 
the  se  :ond  one  is  two-dimensional  and  deals  with  a  partial  differential 
equation.  Regular  perturbations  are  used,  in  both  cases,  to  find  solutions 


far  and  near  singularities  and  the  Rayleigh-Ritz  method  is  used  then  to 
match  the  two  solutions  and  find  the  unknown  constants  they  contain. 
Even  though  both  problems  are  linear,  the  procedure  is  by  no  means  re¬ 
stricted  to  problems  of  this  class. 


1.  DEFLECTION  OF  A  SLIGHTLY  STIFF  MEMBRANE 

Suppose  that  a  pressurized  membrane  (Fig.  1),  instead  of  being  per¬ 
fectly  flexible,  has  a  slight  bending  stiffness.  Without  loss  of  generality 
the  governing  equation  for  the  deflection  w  may  be  assumed  to  have  the 
following  form  £4] : 


£ 


1 


(1) 


The  boundary  conditions  for  clamped  edges  are 


w 


at  x  =  ±  1 


(2) 


The  first  term  in  eq.  (1)  represents  the  effect  of  bending  rigidity  (e 
being  a  small  parameter,  «  1).  The  second  term  describes  the  effect  of 
membrane  forces,  and  the  third  term  the  effect  of  transverse  loading, 
which  is  taken  to  be  uniform  for  simplicity. 

If  one  tries  to  solve  the  problem  by  using  a  regular  perturbation 


scheme,  such  as 


3 


the  system  of  second-order  equations  for  wq  •  wi»  ***  is  obtained,  after 
Introduction  of  eq.  (3)  In  eq.  (1),  and  not  all  four  boundary  conditions 
can  be  satisfied.  This  is  a  typical  problem  encountered  in  so-called 
singular  perturbations. 

It  is  clear  that  the  difficulties  originate  at  the  edges  where  the 
first  term  in  eq.  (1)  is  of  importance.  The  solution  (3)  is  otherwise 
valid  far  away  from  x  ■  ±  1  .  Therefore,  away  from  the  edges  the  so- 
called  global  solution  is 

wG  *=  w0  =  f  d  "  x2)  +  0(e) 

In  order  to  obtain  the  solution  which  is  valid  at  the  edges,  we  fol¬ 
low  the  general  boundary- layer  procedure  [4]  of  shifting  the  origin  to 
the  region  of  nonuniformity  and  then  magnifying  the  scale  by  introducing 
a  new  coordinate  X  instead  of  x  ,  as  follows. 

Due  to  symmetry  we  consider  the  left-hand  edge  only  (at  x  -  -1) 
and  introduce  the  new  coordinate 


X 


x  +  1 


(5) 


where  the  magnification  factor  1/e  was  chosen  in  such  a  way  that  in 
the  new  governing  equation,  given  below,  stiffness  at  the  edge  balances 
tension,  rather  than  pressure  [4],  Therefore,  after  introducing  eq.  (5) 
into  eq.  (1),  the  governing  equation,  valid  near  x  ■  -1  ,  is 


. 4  .2 

d  w  _  d  w 

dX*  dX2 


(6) 
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The  solution  of  this  equation  is  of  the  form 

vL  -  AeX  +  Be"X  +  CX  +  D  (7) 

After  the  exponentially  growing  term  in  eq.  (7)  is  rejected  (on  physical 
grounds)  and  the  boundary  conditions  (at  x  •  -1)  satisfied,  the  locally 
valid  solution  is 

wL  -  C(X  -  1  +  e”X)  +  0(e)  (8) 

The  constant  C  is  unknown  and  is  to  be  found  next.  This  is,  then,  a 
typical  situation  in  which  a  so-called  matched-asymptotic  method  is  used. 
But  we  now  want  to  show  that,  in  order  to  find  C  ,  the  calculus  of  var¬ 
iation  may  be  used  to  good  advantage  also. 

The  idea  is  as  follows.  The  local  solution  is  valid  in  the 

region  -1  s  x  S  x^  ,  where  xA  is  a  certain,  still  unknown,  point  close 

!< 

to  the  edge:  xA  «  -1  +  0(e)  ,  Fig.  1.  The  global  solution  wG  is  valid 
in  the  region  S  *  <  0  .  (We  consider  only  the  left-hand  side  of  the 
membrane.)  The  local  solution  w^  satisfies  the  boundary  conditions  at 
x  -  -1  exactly.  Therefore,  if  one  forms  the  variational  statement  of  the 
problem  and  uses  as  trial  solution  w^  on  -1  <  x  <  xA  and  wG  on 
XA  i  *  -  0  all  the  boundary  conditions  will  be  satisfied  and  the  so-called 
boundary  terms  in  the  variational  statement  will  vanish.  In  that  case 
the  necessary  calculations  reduce  to  a  simple  integration  as  described 
below. 

It  is  well  known  that  the  calculus  of  variations  provides  as  good 
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results  as  the  trial  functions  Is  a  good  approximation  to  the  exact  solu¬ 
tion.  Since  Wg  and  w^  may  be  found  to  any  order  of  accuracy,  one  is 
to  expect  that  the  calculus  of  variation  should  be  able  to  provide  the 
correct  result  for  the  constant  C  . 

The  variational  statement  for  the  system  (1)  and  (2)  is  [5]: 


Here,  U  is  the  potential  energy  of  the  membrane,  6  is  the  variational 
symbol  and  B.T.  stands  for  "boundary  terms".  We  can  write  eq.  (9)  as 


<5  (7}  =6  /  £(w)dx  +  B.T.  =  0 

x  =  -1 


(10) 


where  the  operator  £  is  given  by 


c  a2  2  .  ,  2 
f  s  7(75)  +*(&)  -1 

dx 


2  '  dx  • 


(11) 


The  last  integral,  eq.  (10),  may  be  split  up  into  two  parts  in  terms 
of  local  and  global  solutions  as 

„  XA  0 

<5  ("2)  *  5  /  2(wL)dx  +  /  £ (Wg)dx  -  0  (12] 

X  »  -1  X  ■  X 

A 

Note  that  the  boundary  terms  (B.T.)  vanish,  since  w  is  used  as  the 

U 

trial  solution  at  x  *  -1  .  Note  also  that  there  are  now  two  unknowns  in 
the  problem;  xA  ,  the  place  of  matching  w^  and  wfi  and  the  constant  C 
which  is  hidden  in  w  in  eq.  (12). 

L 


u 

C  ■  e  .  The  same  result  for  C  Is  given  in  [4]. 

Note  that  C  and  x.  can  not  be  more  accurate  than  to  order  of 

A 

u 

e  because  we  started  calculations  with  w_  and  wT  which  are  accurate 

b  L 

up  to  order  smaller  than  e  . 

After  w  and  w,  are  known,  the  composite  solution  valid  every- 
b  L 

where  may  be  found  by  using  Erdely's  scheme,  described  in,  say  [11.  It 
is  given  in  [4]  as: 


w  «  \  (1  -  x2)  -  ^  [l  -  exp  [-  ^ p j] 

2.  TORSION  OF  SHAFT  WITH  KEYWAY 

As  a  second  example  of  the  application  of  the  calculus  of  variations 
in  singular  perturbations  we  consider  a  problem  of  determining  the  stress 
concentration  factor  for  a  circular  cylindrical  shaft  in  torsion  with  a 
semi-circular  keyway  cut  along  its  length,  Fig.  2.  The  exact  solution 
for  this  problem  may  be  found  in  [5]  and  [6],  while  the  same  problem  is 
treated  in  [4]  by  using  matched  asymptotic  expansions. 

Mathematically,  the  problem  reduces  to  that  of  finding  the  stress 
function  <J>  such  that  it  satisfies  the  equation 


2  2 

L±  +  1_±  „  _i 
2  2 
3x  3t 


(15) 


over  the  cross-section  of  the  shaft,  and  the  boundary  condition  d>  *  0 

along  the  boundaries  of  the  cross-section.  Fig.  2.  Then  the  maximum  shear 

3<t> 


stress  t 


3y 


is  found  at  the  root  of  the  keyway  and  compared  to  the 


corresponding  maximum  stress  of  the  circular  cross-section  (without  the 


keyway) . 


We  first  find  the  global  solution  <f>  for  a  circular  cross-section 
(without  the  keyway).  It  obeys  the  differential  equation  (15),  subject  to 
the  boundary  condition  $g(r  =  1)  =  0  .  The  solution  of  this  problem  is 
easily  found  to  be 


(16) 


In  order  to  find  the  local  solution,  valid 
of  the  keyway  (Fig.  3),  we  follow  the  same  idea 
First,  the  magnified  coordinates  are  introduced 


in  the  immediate  vicinity 
as  in  the  previous  example, 
as,  [4]: 


X 


x 

e 


Y 


V  +  1 
e 


(17) 


After  introducing  eq.  (17)  into  eq.  (15)  the  governing  equation  becomes 


2 

-£ 


(18) 


subject  to  the  boundary  condition  <j>  *  0  along  the  boundary  shown  in 

Li 

Fig.  3.  This  local  solution  is 


♦l  -  c(y 


+  0  (e  2 ) 


(19) 


The  problem  reduces  now  to  that  of  finding  the  unknown  constant  C 
in  eq.  (19).  In  order  to  find  C  we  again  make  use  of  the  Rayleigh-Ritz 


procedure  as  in  the  previous  problem. 

Before  performing  the  matching,  we  introduce  polar  global  (r,9) 


9 


and  local  (R,iJ>)  coordinates  (Fig.  2  and  Fig.  3): 


2  2  2 

x  =  r  sin  0  ,  y  =  -r  cos  6,r  =  x  +  y  ,0  =  arc  tang  (-  —  ) 


X  =  R  sin  ip  ,  Y  =  R  cos  \p  ,  R^  =  X^  +  ,  ip  =  arc  tg  ^ 


We  then  formulate  the  variational  statement  for  the  complete  problem  (with 
the  keyway),  and  it  is  given  in  [5]  as 


fiU  =  6 


2$  dx  dy  +  B.T.  =  6  f /  i’ (<)>)dx  dy 

D 


+  B.T.  =  0 


(20) 


where  D  is  the  area  of  interest  and  B.T.  stands,  as  before,  for 
"boundary  terms". 

To  assure  that  4>  satisfies  the  boundary  conditions  everywhere,  we 
write  the  variational  statement  in  the  form 

5U  =  6  //  £(<}>  )dx  dy  +  6  //  i'  (<f>  )dx  dy  E  SU~  +  5U  (21) 

°G  °L 

where  S'  (<J)  is  defined  by  eq .  (20)  and 

*i<v  ■  ( H  )2  +  (14  )2  +  °<t2> 

<p  therefore  satisfies  the  global  boundary  condition  along  the  circle 
G 

(r  =  1)  while  satisfies  the  local  boundary  condition  along  D 

L  L. 


r%*  * 


(Fig.  3). 
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We  now  switch  to  the  polar  coordinate  introduced  previously,  such 

that 


(  —  )2  +  (~  )2  =  (~  )2  +  !;(  —  )2 

3x  3y  3r  r  3$ 

(  A)2+  (  l^)2  +±2  I1-)" 

3X  3Y  3R  R  3t|» 

Also  we  will  slightly  modify  regions  D  and  D  on  which  the  Integra- 

b  L 

tion  in  eq.  (21)  will  actually  be  performed.  Due  to  the  symmetry  of  the 
problem,  we  are  to  integrate  only  across  one-half  of  the  cross-section. 

We  assume  that  the  global  solution  is  valid  inside  the  region  CDEGC 
(Fig.  4),  while  the  local  solution  is  valid  inside  the  region  ABED  .  The 
two  regions  are  divided  by  the  circle  of  radius  k  (of  order  e  and  still 
unknown)  whose  equation  is 

2  2  k 

r  =  cos  0  -  (k  -  sin  8) 

and  for  small  0  is  given  by 


When  determining  U  we  will  first  perform  the  integration  across 
the  semicircle  CFG  and  subtract  the  corresponding  integral  along  DFE  . 
But  the  first  integral  will  be  a  constant,  say  Cj  ,  and  will  not  depend 


on  the  unknowns  (k  and  C) .  Therefore,  the  variation  of  that  integral 


11 


in  Ug  will  be  equal  to  zero.  Thus,  the  integration  necessary  to  find 


Ug  will  only  be  performed  across  the  region  FEDF  ,  where 
0  <  6  s  k  =  tg0  and  r  changes  from  r  *  1  -  (k  -6)’  to  r 
Therefore,  we  find  that 

2 


1  . 


0°  ’  '  ABED  H  »')  +  M  >>)  "  2*4 


where  the  Jacobian  J  is  given  by 

u 


JG  dx  dy  + 


(23) 


JG  = 


3x 

3x 

3r 

30 

ix 

lx 

3r 

30 

(24) 


After  the  integration  with  respect  to  r  is  performed  in  eq.  (23)  and  some 
higher  order  terms  disregarded,  the  following  result  is  obtained 


K 


-  j  f  (k2  -  e2)h  de 


-  V2  +  r 
8  *  L1 


(25) 


The  same  procedure  is  now  repeated  with  the  local  solution  valid  in- 

7T  k. 

side  the  region  ABED  where  0  a  ip  <  j  and  1  i  R  <  —  .  Therefore,  one 
has  to  find  U  from 

L 


\2 


d4> 


(26) 


ABED 


where  the  Jacobian  J^  =  -R  may  be  found,  by  using  local  coordinates  and 
replacing  them  in  eq.  (24)  instead  of  global  ones.  is  given  by  eq. 

(19)  and  in  polar  coordinates  as: 


12 


C(R  -  -)  cos 


(27) 


After  the  necessary  integration  in  (26)  is  completed,  one  finds  that 


L 


(28) 


Knowing  U„  and  U.  from  eq.  (25)  and  (28),  the  total  potential 
b  L 

energy  may  be  found  as 


U(k,C)  =  UL  +  UG 


(29) 


Performing  the  variation  of  U  with  respect  to  two  unknowns  (k  and  C) , 
one  obtains  the  following  equations 


One  can  again  use  the  perturbation  scheme  such  as  in  eq.  (14)  to  find 
2  12 

that  k  =  e  +  0(e  )  and  C*  j£+0(e  )  .  The  same  result  for  C  was 
obtained  in  [4]  and  that  is  in  essence  what  we  wanted  to  show. 

The  composite  expansion  is  found  again  by  using  Erdely's  scheme  [1] 
and  is  given  by 


2 


x 


yf  i.  - 

+  (1  +  V)2 
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The  maximum  shearing  stress  (at  the  bottom  of  the  keyway)  is  given  as 


kK  9y 


t„  *  t-1 2 3 4 5 6:  *  1  +  0(e) 


while  in  the  case  when  there  is  no  keyway  the  corresponding  stress  is 
3*G  1 

t  ,  =  — —  =  7t  .  Therefore,  the  stress  concentration  factor  is 
wk  3y  2 

Tk/Twk  =  2- 
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FIGURE  CAPTIONS 


Fig.  1  Pressurized  Membrane 

Fig.  2  Cylindrical  Shaft  in 

Fig.  3  Vicinity  of  Keyway 

Fig.  4  Part  of  Shaft 
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